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ABSTRACT 

We present the AGN, star-forming, and morphological properties of a sample of 13 MIR-luminous (/24 > 
700 /iJy) IR-bright/optically-faint galaxies (IRBGs, /24//R > 1000). While these z ^ 2 sources were drawn 
from deep Chandra fields with > 200 ks X-ray coverage, only 7 are formally detected in the X-ray and four 
lack X-ray emission at even the 2(t level. Spitzer IRS spectra, however, confirm that all of the sources are 
AGN-dominated in the mid-IR, although half have detectable PAH emission responsible for ^25% of their 
mid-infrared flux density. When combined with other samples, this indicates that at least 30-40% of luminous 
IRBGs have star-formation rates in the ULIRG range 100-2000Mq yr~'). X-ray hardness ratios and 
MIR to X-ray luminosity ratios indicate that all members of the sample contain heavily X-ray obscured AGN, 
80% of which are candidates to be Compton-thick. Furthermore, the mean X-ray luminosity of the sample, 
log L2-iokev(ergs s"')^ 44.6, indicates that these IRBGs are Type 2 QSOs, at least from the X-ray perspective. 
While those sources most heavily obscured in the X-ray are also those most likely to display strong silicate 
absorption in the mid-IR, silicate absorption does not always accompany X-ray obscuration. Finally, ^ 70% 
of the IRBGs are merger candidates, a rate consistent with that of sub-mm galaxies (SMGs), although SMGs 
appear to be physically larger than IRBGs. These characteristics are consistent with the proposal that these 
objects represent a later, AGN-dominated, and more relaxed evolutionary stage following soon after the star- 
formation-dominated one represented by the SMGs. 

Subject headings: galaxies: active — infrared: galaxies — X-rays: galaxies 



1. INTRODUCTION 

Studies of heavily-obscured star-forming galaxies and ac- 
tive galactic nuclei (AGN) in the distant Universe have lagged 
behind those of their unobscured counterparts, due largely to 
their being extremely faint in the optical and UV. This op- 
tically faintness, however, need not be a limitation and can 
instead be used as a selection criterion, as dust-enshrouded 
sources faint at short wavelengths should be comparably 
bright in the infrared where their absorbed radiation is re- 
emitted. One might therefore expect sources with bright in- 
frared emission yet faint optical magnitudes (/24//R > 1000) 
to be ideal obscured galaxy/AGN candidates, and indeed, 
such s ources comprised so me of the first Spitzer Space Tele- 
scope (IWerner et al.ll200l targets. 

Initial studies of the brightest (/24 > 700 /iJy) IR- 
bright/optically-faint galaxies (referred to here as IRBGs) tar- 
geted for Spitzer IRS mid-infrared spectroscopic ( Houck et"an 
120041) follow-up observations indicated that most lie at z ~ 2 
and have either featureless spectra or spectra dominated by 
silicate absorption^, properties indicative of AGN activity (e.g., 
iHoucketalJ I200I iYan etalJ llOOSt iWeedman et all 12006a). 
Their optical faintness was therefore attributed both to dis- 
tance as well as to obscuration by the dust surrounding an 
AGN's central engine, although more recent studies suggest 
that the AGN's host galaxy may also contribute significantly 

' Space Telescope Science Institute, 3700 San Martin Diive, Baltimore, 
MD 21218; Giacconi Fellow; donley@stsci.edu 

^ Steward Observatory, University of Arizona, 933 North Cherry Av- 
enue, Tucson, AZ 85721 

^ Department of Physics, Durham University, Durham DHl 3LE, UK 
Departamento de Astrofi'sica y CC. de la Atmosfera, Facultad de CC. 
Fi'sicas, Universidad Complutense de Madrid, 28040 Madrid, Spain 

' Associate Astronomer at Steward Observatory, The University of Ari- 



to the observed extinction (iBrand et al.l l2007t iPolletta et alJ 

^oi. 

Subsequent multiwavelength studies have since confirmed 
the AGN nature of the most luminous IRBGs, although 
they reach differing conclusions concerning the relative 
importance of star-formation and AGN activity amongst 
fainter 0/24 < 700 ii]y) IRBGs (e.g.. iDev et al] mm 
Donlev et al. 2008; Georgantopoulos et al. 2008; Fioreetal 
2008; ' Pope et aL ,2008a:" Desai et al.. .2009; Treister et al 



2009; Fiore et aL feoogj) . Furthermore, they have shown that 
the space density of the most luminous IRBGs is compara- 
ble to that of luminous, unobscured AGN at z = 2 jDev et alj 
[20081) . If these sources are in fact heavily obscured AGN, 
they therefore represent an important phase in the growth of 
supermassive black holes during the era in which both star- 
formation and AGN activity peaked. 

To constrain the level of obscuration in these luminous 
AGN, studies turned to the X-ray. Unfortunately, the large- 
area survey fields from which the luminous, and therefore 
rare, IRBGs were initially selected (e.g. the Spitzer First- 
Look Survey, the NOAO Deep Wide Field Survey, and the 
SWIRE Survey) and then foUowed-up with IRS have minimal 
X-ray c overage capable of detecting only the brightest AGN 
jDev e t al. 2008). Nonetheless, studies in these fields suggest 
that at least ^ 50% and perhaps as many as 95% of X-ray- 
detect ed IRBGs with /24 > 1 .3 mJy are X-ray-obscured AGN 

S Lanzuisi et a Jl2 009). so me of which may b e Compton-thick 
Polletta et al 



I2OO6: Alexander et al.ll2M8bl) . 



Alternatively, other studies have focused on samples of IR- 
BGs selected in fields with deeper X-ray coverage, such as 
the Chandra Deep Fields and COSMOS. The deepest X- 
ray fields, however, are also the smallest, so those studies 
most suited to constrain X-ray obscuration are also those in 
which the focus generally shifts to the fainter and more nu- 
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TABLE 1 
Observations 



ID 


Field 


Program 


Observation 


Ramp Duration (s) 


Number of Cycles 


Total 






ID " 


Date 


LLl'' 


LL2'^ 


LLl'' 


LL2'^ 


Integration (s)"* 








zuuo iepi ij 


120 


120 


10 


10 




IRBG2 


CDF-S 


30419 


2007 Aug 31 


120 


120 


23 


23 


11040 


IRBG3 


E-CDFS 


30419 


2006 Sept 14 


120 


120 


27 


27 


12960 


IRBG4 


CDF-S 


30419 


2007 Aug 31 


120 


120 


15 


15 


7200 


IRBG5 


CDF-N 


20456= 


2006 Apr 24 


120 


120 


6 


6 


2880 


IRBG6 


CDF-N 


20733 


2006 May 26 


30 


30 


80 


80 


12880 


IRBG7 


CDF-N 


20456' 


2006 Apr 22 


120 


120 


6 


6 


9600 


IRBG8 


EGS 


30419 


2007 Apr 27 


120 


120 


12 


12 


5760 


IRBG9 


EGS 


30419 


2006 Jun 26 


120 


120 


4 


4 


1920 


IRBGIO 


EGS 


30419 


2007 Apr 27 


120 


120 


7 


7 


3360 


IRBGII 


EGS 


30419 


2007 Mar 25 


120 


120 


28 


28 


13440 


IRBG12 


EGS 


30419 


2007 Mar 25 


120 


120 


18 


18 


8640 


IRBG13 


EGS 


30419 


2007 Apr 27 


120 


120 


6 


6 


2880 



" The Pi's of programs 30419, 20456, and 20733 are G. H. Rieke, R. Chary, and C. M. Urry, respectively. 
Long Low 1st order, 19.5-38.0 
Long Low 2nd order, 14.0-21.3 /im 

The IRS nodding peiformed in Standard Staring mode results in two spectra per cycl e, each with an exposure time equal to the ramp duration. 
The IRS spectrum for this source was previously published in lPope et alj )2008bl) and lMurphv et alj (200% 
' The IRS spectrum for this source was previously published in lMurphv et alj 120091) 



merous population of IRBGs whose nature rem ains contro- 
versial. In the larger COSMOS field, however, iFiore et al] 
([2009) find that 40% of IRBGs with moderate flux densi- 
ties of /24 > 550 /iJy are detected in the X-ray with a mean 
rest-frame obscured luminosity of log L2-iokev(ergs s"')= 43.5 
and an X-ray hardness ratio (HR=(H-S)/(Hh-S), where H=1.5- 
6 keV and S=0.3-1.5 keV) of 0.50, consistent with obscured 
yet Compton-thin absorption. The stacked X-ray signal from 
the remaining 60% of the sample has a similar HR of 0.53, 
although its interpretation is dependent on a number of fac- 
tors including the assumed underlying column density distri- 
bution, the intrinsic MIR/X-ray luminosity ratio, and most 
importantly, the assumed photon index of emission from 
star formation. Assuming a star-forming X-ray photon in- 
dex (r) of 1.9, IFiore et alj ( 120091) conclude that 94% of the 
X-ray non-detected, f24 > 550 /iJy IRBGs are heavily ob- 
scured, and likely Compton-thick, AGN. If a harder X-ray 
photon index of F ~ 1 .0 - 1 .4 (as observed in th e starburst- 
dominated ULIRGS of 'Franc eschini eta n (l2003h . |Ptak et alj 
([2003), and lTeng et a l. ( 2005)) is assumed, however, the ob- 
scured AGN fraction of this moderate-luminosity sample falls 
significantly. 

To bridge the gap between these infrared and X-ray studies 
and therefore better constrain the luminosity, obscuration, and 
power sources of these cosmologically interesting AGN, we 
have obtained IRS spectra of a sample of IRBGs selected in 
deep X-ray fields with effective exposures of Tx > 200 ks. Be- 
cause we do not require that our targets be detected in the X- 
ray, however, we do not bias ou r sample tow ards the brightest 
or least-obscured galaxies (e.g.. lBrand et alj |2008). The sam- 
ple discussed below is therefore the first uniformly-selected 
sample of IRBGs with both deep X-ray and IRS coverage. 
Because they were selected from deep multiwavelength fields, 
many of the sources also have Hubble Space Telescope (HST) 
ACS imaging, enabling a morphological study of this unbi- 
ased sample and a comparison to sub-mm galaxies (SMGs), a 
population of z 2 star-forming galaxies proposed by some 
to be the merger-induced progenitors of luminous IRBGs 
(e.g.,|De v et al. 2008; Pope et al. 2008a; Coppin et al. 2010.; 
iNaravanan et alJI2009bllal;'lBrodwin et alJI2008l) . 



The paper is organized as follows. In §2, we discuss the 
sample selection, observations, and data reduction. The op- 
tical and infrared photometric properties of the sources are 
then examined in §3. In §4, we present the IRS spectra and 
discuss the AGN and star-forming contribution to the MIR 
emission of these sources. The X-ray emission is discussed in 
§5, as is the agreement between the X-ray and IR properties, 
and the star-formation rates of the sources are discussed in 
§6. In §7 we summarize the morphological properties of the 
IRBGs' hosts and compare them to those of the SMGs. The 
discussion follows in §8, and we then summarize our con- 
clusions in §9. Throughout the paper, we assume the follow- 
ing cosmology:(f}„,17A,//o)=(0.27, 0.73, 70.5 km s"' Mpc"'), 
and quote all magnitudes in the AB system unless otherwise 
noted. 

2. SAMPLE SELECTION AND OBSERVATIONS 

To ensure deep X-ray coverage, we selected the IRBG sam- 
ple from the Chandra Deep Fields North and South (CDF-N, 
CDF-S, T; - 2 Ms; Alexander et al. 2003; Luo et al. 2008), 
the Ext ended Chandra D eep Field South (E-CDFS, 7^ ^ 
250 ks; lLehmeretn]|2005i). and the Extended Groth Strip 
(EGS, n ^ 200 ks; iLaird et alj|2009h . flie combined area of 
which is 0.96 deg^. We chose for IRS follow-up sources with 
extreme IR/optical flux ratios (/24 //r > 1000) typical of those 
in Houck et dj (120051) . lYan et alj (12005). and iWeedman etal] 
f2006 a.b) based on the Spitzer MIPS 24 /xm and optical 
photometry compiled in lPerez-Gonzalez et al.l (120081) and the 
UCM Extragalactic Databas^. In addition, we required flux 
densities in excess of /24 > 700 /iJy to guarantee that IRS 
spectra could be obtained in a reasonable amount of time. The 
final sample is comprised of 13 IRBGs. 

The 3 sources in the CDF-N were observed with IRS as 
part of programs 20456 (PI Chary) and 20733 (PI Urry), 
an d the spectra of tw o (I RBG5 and IRBG7) c an be found 
in Pop e et alj ( l2008bl) and Mu rphy etalj (12009'). For the re- 
maining 10 sources, we obtained 1st and 2nd order spectra 
with the IRS Long-Low (LL) module (program 30419, PI 

* http://guaix.fis.ucm.es/~pgperez/Proyectos/ucmcsdatabase.en.htmI 
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TABLE 2 

IR-Bright/Optically-Faint Sample 



ID RA(J2000) DEC(J2000) z^pec" ziRs'' /24 R(AB) /24//R "iRAc"^ QlRs"^ AGN"*" /24.sf'= rgj^m' U^im" 













yh^-'y ) 


















IRBGl 


03:31:46.6 


-27:45:53.0 






1069 


25.3 


3961 


-1.92 


-2.08 


96 ±2 


41 


0.0 


45.6 


IRBG2 


03:31:58.3 


-27:50:42.2 




2.28± 0.06 


851 


25.2 


2934 


-1.86 


-2.43 


94 ±5 


99 


0.5 


45.6 


IRBG3 


03:32:10.5 


-28:01:09.2 




1.62± 0.02 


751 


25.1 


2208 


-1.46 


-2.06 


71 ±4 


286 


0.4 


45.0 


IRBG4 


03:32:37.8 


-27:52:12.5 


1.60 




1066 


24.8 


2370 


-1.78 


-2.07 


100 ± 3 





0.0 


45.2 


IRBG5 


12:36:00.2 


62:10:47.3 


2.00 


2.00± 0.03 


1144 


24.0 


1204 


-2.13 


-2.06 


84 ±3 


293 


0.1 


45.4 


IRBG6 


12:36:35.6 


62:14:23.9 


2.02 


2.04±0.004 


1497 


23.8 


1426 


-2.00 


-2.04 


85 ± 1 


351 


0.0 


45.6 


IRBG7 


12:37:26.5 


62:20:27.0 




1.76± 0.02 


883 


25.0 


2324 


-2.21 


-1.81 


> 72*' 


272 


5.2 


45.4 


IRBG8 


14:15:38.1 


52:18:52.9 






1217 


24.9 


2967 




-1.29 


90 ±6 


122 


0.1 


45.7 


IRBG9 


14:16:35.4 


52:12:35.2 






2299 


23.7 


1999 


-2.53 


-1.92 


92 ±4 


327 


0.2 


45.9 


IRBGIO 


14:18:51.5 


52:48:35.9 


1.92 


1.99± 0.02 


1686 


24.0 


1812 




-2.32 


88 ± 3 


394 


0.3 


45.6 


IRBGl 1 


14:20:08.6 


53:01:10.1 




2.07± 0.02 


701 


24.7 


1491 




-1.94 


92 ±9 


57 


5.0 


45.5 


IRBG12 


14:20:13.0 


52:55:33.3 




1.04± 0.01 


948 


24.8 


2106 




-2.64 


78 ± 3 


218 


0.1 


44.3 


IRBG13 


14:20:21.9 


52:55:12.1 




2.60± 0.06 


1786 


24.1 


2152 




-0.58 


79 ±5 


440 


1.4 


46.1 



References. — Spectroscopic redshifts from lSzokoIv et alj f200l IRBG41 ISwinbank et all f20Ql IRBG5, IRBG6), and lDavis etaH <2007l IRBGIO) 
" Optical/NIR spectroscopic redshift from the literature 
IRS-derived redshift (this work); see §4 

IRAC (3.6-8.0 fim) and IRS (14-37 ^m) power-law slopes, where oc i/" 
'I AGN contribution to the MIR (14-37 ^m) flux density 

Star-forming contribution to the observed 24 /im flux, as determined by the spectral fit 
' Optical depth of best-fit AGN continuum at 9.7/^m (see §4) 
^ Rest-frame 6 fim luminosity of the AGN compononent (see §5.2) 

IRS spectrum not fit by default model. See §4.1 for details. 



TABLE 3 
IR Selection Criteria 



ID 


H05 


Y05 


W06 


D08 


F08 


G08 


P08 L09 


PL 


IRBGl 


X 




X 


X 


X 


X 




X 


IRBG2 


X 






X 


X 


X 


X 


X 


IRBG3 


X 






X 




X 




X 


IRBG4 


X 




X 


X 


X 


X 


X 


X 


IRBG5 




X 




X 






X 


X 


IRBG6 








X 


X 


X 


X 


X 


IRBG7 


X 






X 




X 


X 


X 


IRBG8 


X 




X 


X 


X 


X 


X 




IRBG9 








X 


X 


X 


X 


X 


IRBGIO 




X 




X 










IRBGl 1 




X 




X 










IRBGl 2 


X 


X 




X 


X 


X 






IRBGl 3 






X 


X 






X X 





ual background, rogue pixels (pixels with high dark current 
and/or photon responsivity), and cosmic rays, and takes into 
account all frames produced by the IRS/SSC pipeline. The 
background and noise images are produced by masking the 
target spectrum (and any serendipitous spectra) on each frame 
and then coadding the resulting frames. A biweight statistical 
estimator is then used iteratively and interactively to minimize 
the contamination from deviant pixels, and pixels that deviate 
by more than Scr from the mean local value are flagged as 
rogue. After rejecting any cosmic ray events and taking the 
spectral distortion into account, the spectra are then weighted 
by the PSF, optimally extracted, and coadded to produce the 
final spectrum. 



References —H05 : Houck et al. 1 200^, Y05: lYan et al] f200l. W06: 
IWeedman et al.' (2006a"b), DOS: Dev et al. (200^ (DOGS), F08: Fiore et alj 
12008), G08: Georgantopoulos et al. (2008), P08: PoUetta et al. (20oC 
L09:.Lanzuisi et al..t2009.) ( EDOGS). PL (power-law): .Alonso-Herrero et al] 
<2006t) : IDonlev et alj <200l) 

Note. — An 'x' indicates that a source would be selected via the listed 
criteria and a '— ' indicates that a source lacked the neccessary data to deter- 
mine if it would be selected. 

Rieke). Observation details are given in Table 1. The re- 
sulting wavelength range of 14-37 pim guarantees coverage 
of the 9.7 /im silicate absorption feature at z = 0.4-2.8 and 
the 7.7 fim aromatic (hereafter poly cyclic aromatic hydrocar- 
bon (PAH)) feature at z = 0.8-3.8. Of the 13 sources in our 
sample, 4 have optical/NlR spectroscop ic redshifts from the 
hterature, placing the m at 1 .6 < z < 2.0 (ISzokoly et al.ll2004t 
ISwinbanket ani2004tlDavis et al.ll2007l) . 

2.1. Data Reduction 

The IRS data were reduced using the IrsLow package de- 
veloped by D. Fadda to accurately measu re low-resolution 
spectra of faint, high redshift sources (see F adda et al.ll20To[ 
for more details). Briefly, this package corrects for resid- 



3. OPTIC AL/IR PHOTOMETRIC PROPERTIES 

We list in Table 2 the optical and IR characteristics of the 
IRBGs. The 13 sources in our sample have 24 /im flux den- 
sities ranging from 701 to 2299 fiiy (median /24 = 1069 fiJy), 
R-band magnitudes ranging from 23.7 to 25.3 AB (median 
R = 24.8 AB), and 24 fim to R-band flux ratios ranging from 
1204 to 3961 (median /24//R = 2152). 

To place these sources in the context of recent stud- 
ies, we also show in Table 3 the IR-excess selection 
criteria met by each of the 13 IRBGs. A ll sources 
meet the DOG (dust-obscured galaxy) criteria of iDey et alJ 
(I200I f24/fR > 1000 /24 > 0.3rn J y), an d aU but 2 meet 
one or more of the iHouck et al.l (120051 /24 > 0.75inJy, 
iJvega > 24.5'). lYan et al.l(l2005L los(v fJ24u m) /lyf, AS um)) > 
0.5, log(t / f,A24-nm)/iy fJO.lum)) > 1.0), or ' We edman etaP 
(I2006allbl /24 > l.OmJy, T^vega > 23.9) selection crite- 
ria. Fu rthermore, n ine so urces meet the IR-excess cri- 
teria of jRore et al.' ('20 01 hAlfR > 1000,(7? - /:)vega > 
4.5) or [Georgantopoulos etal] (|20M /24//R > 1000,;;ab- 
wJs.e/im.AB > 3.7), and eight meet the combined MIPS/IRAC 
criteria defined bv .Polletta et al.. (.2008.) to select highly lumi- 



4 



DONLEY ET AL 



\est(Mm) 



1.0 10.0 



.IRBG1 




..f 


P.I. 


1.0 


10.0 


0.1 1.0 


10.0 


IRBG6 

z = 2.02 

•• 














Pi: 


1.0 


10.0 


0.1 1.0 


10.0 


ikBGTf 

z = 2.07 




■ 





.ifRBG'2 ■' 

z = 2.28 






Pi 



.0 10.0 
0.1 1.0 10.0 



1.0 



10.0 






0.1 1.0 


10.0 




15 




.lRBG'3"' 




15 






z=1.62 






20 


o 




* 




25 


o 


/ 
r 




25 






, . ..P.I 


30 






1.0 


1 0.0 




o 








15 




IRBG8 


^ • 


1 5 


20 




• • 




20 












25 




• 

• 




25 




o 


f 










1.0 


10.0 








0.1 1. 


10.0 




15 


o 


'lftBGi'3 




15 






z = 2.60 






20 


o 


• • 




20 


25 








25 




o 












1.0 


10.0 





1.0 10.0 




1.0 10.0 



1.0 10.0 



■IRBG9 








m 


R,L 



1.0 10.0 



.IRBG5 ' 

z = 2.00 




. * 


Ri 


1.0 


10.0 


0.1 1.0 


10.0 


TRBGirci 

.z=1.92 


X • 






1.0 


10.0 



15 




20 




25 






0) 




"O 


15 


D 




C 


20 


cn 

D 




E 


25 






CD 




< 



Fig. 1. — Broad-band SEDs of the 13 IRBGs. The majority show power-law SEDs that extend from the near to mid-IR. Those sources that meet the formal 
IRAC power-law criterion (see §3.1) are indicated by a 'pi' in the lower right comer. The best fits to the IRS spectra (see §4) are shown in solid (red) Unes and 
the template of NGC 6090 is over-plotted (in green) on the SED of IRBG12 to confirm its low redshift (z ~ 1; see §4.1). 



nous and heavily obscured AGN. 



3.1. SEDs 



The broad-band (0.4-70^m) SEDs of the IRBGs are shown 
in Figure 1. We constructed the SEDs by first compil- 
ing apertur e-matched photometric catal ogs using the CDF-S 
data of Marzke et alj a999, j;/z). ,Vandame et al.l (i2001. JK) 
ArnoutsetalJ (I2002L U U„BVRI\ COMBO 17 jWolfetal 



iGiavahsco etal] (12004 bvizJHK), iLe Fevre et a" 



(12004 I), and C rALEX (FUV,N UV), the CDF-N da ta of 
IGiavahsco et"an (|200 4'. bviz) and ICapak et al.1 (12004', UB- 
VRIz'HK'), and the EGS data (UBgVRIzJK) of Villar et aD 
(|2008) and Barro et al. (2010), in preparation. For more de- 
tails on the available datasets and the aperture -matching pro- 
cedure, see lPerez-Gonzalez et all (I2005L l2008l) and the UCM 
Extragalactic Database (see footnote 6). 

While a handful of the source SEDs shown in Figure 1 dis- 
play a weak stellar bump (the broad 1 .6 /im feature that domi- 
nates the SEDs of star-forming galaxies and which tends to be 
correlated with t he presence of PAH emission in IRBGs (e.g. 
iDesai et aTl 12009 )). the majority of the sources in our sam- 
ple have power-law SEDs indicative of AGN activity. When 
we apply the IRAC (3.6-8.0 i-im) power-law selection cri- 
teria used by Alonso-Herrero et al. (20(3^ and iDonley et alJ 
(j2007) to identify AGN-dominated sources, we find that 7 of 
the 13 sources are indeed power-law AGN and one additional 
source (IRBGS) meets the criteria if we lower the required 
chi-squared probability from Px > 0.1 to Px > 0.01 (e.g., 
[Alonso-Herrero et al. 20060. Because the IRAC photome- 
try of IRBGS is questionable due to its position near the edge 
of the IRAC field, the total power-law fraction is ~ 8/12, or 
^ 70%. Our results are therefore in agreement with those 



For consistency with the IRAC power-law selection in 

we assume an 



For consisten cy wim t ne IKAL power-la 

lAlonso-Herrero et al.l 12003) and IDonley et ai] <200% , 
IRAC flux calibration uncertainty of 10% 



of iDey et alJ (l2008h . who find that while only - 20% of IR- 
BGs at /24 = 300 ^Jy have power-law SEDs, the power-law 
fraction rises to ~ 55% at /24 = 700 /iJy (our flux cut) and 
to 70-85% at /24 > 1000 /iJy (the median flux density of 
our sample). Furthermore, all of the power-law sources in our 
sample are extremely red, with a mean slope of a = -1.99 
(where oc i/"), and are the refore strong AGN c andidates 
regardless of redshift (see Do nley et alJl2007ll2008l) . 

3.2. Infra red colors 

A number of infrared color-criteria have been used in the 
literature to discriminate bet ween AGN an d star-formation 
activity (for a review, see IDonley et alJ [2008). While 
the reliability and completeness of pure color/colo r cuts 
is heav ily dependent on a source's redshift. Pop e et aP 
(I2008al) find that for IRBGs (which typically lie at"z - 2), 
/8.o/j,m//4.5/im = 2.0 provides a convenient dividing line be- 
tween star-formation-dominated (/8.o/im/./4.5/im < 2.0) and 
AGN-dominated (./8.o/im//4.5/im > 2.0) emission (also see 
Coppin et al. 2010, who suggest an even lower cut of 1.65). 
As a confirmation of this method, we plot in Figure 2 the ex- 
pected colors of a number of AGN and star-forming galax- 
ies. We find that only AGN fall redward of the cut, as ex- 
pected, although we note that lower-luminosit y AGN (e.g., 
Seyfert galaxies) as well as self-absorbed AGN dTreister et alj 
[2009) could display bluer colors more typical of star-forming 
galaxies. As shown in Figure 2, however, all of our IRBGs 
lie in the red AGN-dominated regime, with 3 sources show- 
ing a significantly higher fiA^ml .h.o^m ratio than the rest. 
iPope et al] ( l2008al) attribute this excess of 24 /im emission ei- 
ther to intrinsic obscuration indicative of AGN activity, which 
enhances the observed-frame 24 /im flux and/or obscures the 
observed-frame 8 /im flux, or to the 7.7 /im PAH feature pass- 
ing into the MIPS 24/im band. We explore these possibilities 
further in §4.3. 

4. IR SPECTRA 

The IRS spectra of our 13 IRBGs are shown in Figure 3. To 
measure redshifts and determine the source of the emission 
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Fig. 2. — MIR colors of the IRBGs, where the numbers give the source ID. 
Overp lotted are the redshifted (z = 1 - 3) templates o f AGN fro mlPoll etta et alj 
(2008', TQSO, Seyfert 2, Seyfert 1.8, Mrk231) and'Assef et al. ( 2008, AGN 
and AGN2) and of purely star-forming ULIRGS from Rieke et al. (2009). 
The two AGN templates with colors blueward of 58.o/54.5 = 2 are the Seyfert 
1.8 and Seyfert 2 templates of Polletta et al. 1 2008). Triangles give the tem- 
plate colors at z = 1 and circles given the colors at z = 3. All of the IRBGs have 
5g,o /54 5 > 2 as expected only for AGN-dominated objects (e.g.. Pope et al. 
l2508ai) . IRBG13, with its rapidly rising IRAC continuum, falls off the plot at 
values of Sg o/S^ s = 12.8 and 524 /58.0 = 8.1. 

(e.g., AGN or star-formation), we simultaneously fit a power- 
law continuum and a star-forming template to each spectrum. 
For the star-forming templates, we adopt the starburst, lu- 
minous infrared galaxy (LIRG), and ultraluminous infrared 
galaxy (ULIRG) templates of Rieke et al. ( 2009). While we 
allow the power-law component to be reddened according to 
the extinction law of lDraind (12003b . in no case do we require 
additional reddening to be applied to the star-forming tem- 
plates. 

Using a chi-squared minimization routine, we then fit the 
spectrum while leaving the template normalizations, power- 
law index, power-law extinction, and redshift as free parame- 
ters. To protect against local minima, we varied the assigned 
starting values for the parameters, and chose the input val- 
ues that produced the fit with the lowest reduced chi-squared 
(if more than one minimum was found). Finally, to place er- 
rors on the fit parameters, we created 1000 simulated spectra 
for each source (based on the measured spectra and their as- 
sociated Gaussian errors) and re-ran our fitting program on 
each. The errors quoted below for the redshifts (§4.2) and 
AGN fractions (§4.3) are taken to be the standard deviations 
in the resulting distributions of fit values. This procedure pro- 
duces reasonable spectral fits for all sources except 1RBG7, 
discussed below, although for 4 objects the lack of spectral 
features prevented direct redshift constraints. 

4.1. IRBG7 

IRBG7 has two spectral features and/or artifacts that impact 
the fit. The first is a sharp emission feature at Aobs ^21 fim, 
or Aemit ^ 7.6^m assuming ou r best-fit redshift of z = 1.76 
(for comparison, iMurphv et al.l i200 9 measure a redshift of 
z = 1.75). This feature, which can be seen in the raw spec- 



tral data, could be strong 7.65 /im [NeVI] emission that is 
partially blended with the 7.7 /im PAH emission. 

The second strange feature in this spectrum is the appar- 
ently blue continuum at A < 1 8 /im, present in both our reduc- 
tion of the IRS data as well as that of Pope et al. (2008b, A. 
Pope 2009, private communication) and Murphy et al. (2009j ). 
If the AGN continuum fit is free to vary over all spectral in- 
dices, it therefore settles on a blue slope of a = 1.69 (where 
fi, cx J/"), despite its extremely red (a = -2.21) IRAC power- 
law continuum. Such behavior, however, appears to be un- 
physical, for while the IRS slope need not match the IRAC 
slope, the two are generally close. In fact, the offset between 
the IRAC spectral slope (aiRAc) and the IRS spectral slope 
(ajRs) for the 7 additional power-law galaxies in our sample 
ranges from 2% to 42%, with a mean (median) value of only 
18% (16%). 

We therefore adopt a revised method to fit the spectrum 
of IRBG7. First, we place an upper limit on the IRS AGN 
continuum slope by adjusting the measured IRAC slope of 
a = -2.21 upwards (bluewards) by 18%, the mean offset be- 
tween the IRAC power-law and IRS slopes. With this as- 
sumed AGN continuum of a = -1.81, we then perform an 
AGN-only fit at A > 18 /im to determine the AGN extinction 
that best fits the profile of the 9.7 /im silicate absorption fea- 
ture. We then set the normalization of the AGN continuum by 
extrapolating the IRAC continuum fit from 8 to 14 //m. When 
we do the same for the other power-law galaxies in the sam- 
ple, we find that the mean offset between the predicted and 
observed 14 /im flux density is 35%. To place a lower limit 
on the AGN contribution, we therefore scale the AGN contri- 
bution downwards until the 14 /im AGN flux density is 35% 
lower than predicted by the IRAC extrapolation, while allow- 
ing the star-forming contribution to vary so as to give the low- 
est reduced chi-squared. The resulting fit under-predicts both 
the short- wavelength emission and the 21 /im emission peak, 
as expected, but now places a robust lower-limit on the AGN's 
contribution to the MIR emission. 

4.2. Redshifts 

Of the 13 sources in our sample, 4 ha ve optical/NIR spec- 
troscop ic redshifts fr om the li terature ( Szokoly et al.l 120041 : 
ISwinbank et al.. 2004i: lDavis et al . 2007). Using the IRS spec- 
tra, we determine reliable redshifts for an additional 6 sources, 
for an overall redshift completeness of ^ 80%. The IRS red- 
shifts derived for 3 of the 4 sources with prior optical/NIR 
redshifts agree to within Az < 0.1 in all cases (The fourth 
source, IRBG4, has a power-law IRS spectrum from which 
no redshift estimate can be made). 

Of the 3 remaining sources, two (IRBGl and IRBG9) have 
power-law spectra and SEDs, and one (IRBG8) has both a 
spectrum dominated by power-law emission as well as ques- 
tionable IRAC photometry due to its position near the edge 
of the IRAC field, thus preventing accurate spectroscopic or 
photometric redshift determination. 

The mean and median redshifts of our sample, z= 1 .89 and 
z = 2.00, respectively, agree quite well with the typic al red- 
shifts o f si milar samples. For in stance, both Weedma n et alj 
(12006a) and 'Houck et al.1 (12005?) find median redshifts of z = 
2.1, Yan et al. (2005) finds a median redshift of z = 2.3, and 
Polletta et al. (2008) finds a median redshift of z = 2.2. 

Finally, with the exception of IRBGl 2, all of the sources 
lie at z > 1.60. The redshift of IRBG12, z = 1.04, is 
therefore quite low in comparison both to our and other 
(e.g., , Weedman et al.. ,2006a; .Houck et al.. .2005; lYan et alj 
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Fig. 3. — IRS spectra of the IRBGs. The green dashed Hne shows the observed contribution from the AGN whereas the blue dotted line shows the observed 
contribution from star-formation. The overall fit is given by the red solid line. 



120051; iPoUetta et al.ll2008"l) samples of luminous IRBGs. This 
source's low redshift, however, appears robust. The reduced 
chi-squared fit has 3 minima, one at z = 0.11, one at z = 1.04, 
and one at z = 3.19. As shown in Figure 1, however, this 
source has a weak 1.6 fim stellar bump that is well-fit at 
optical-NIR rest-frame wavelengths by the z ^ 1 template of 
the starburst galaxy NGC 6090 from Polletta et al. (2006). We 
can therefore rule out the low and high redshift solutions from 
the IRS spectroscopy on the basis of the broad-band SED. 

4.3. Source of power: AGN activity or star-formation? 

The fraction of the overall IRS (14-37 /im) dust-obscured 
flux density arising from the AGN power-law component 
(e.g., the AGN fraction), the best-fit power-law indices, and 
the optical depths at 9.7 /im (for the power-law component) 
are shown in Table 2. On the basis of their MIR (14-37 /im) 
emission, all of the IRBGs in our sample are AGN-dominated, 
with AGN fractions ranging from 71% to 100%. Six of the 
IRBGs are nearly pure power-laws with AGN fractions ex- 
ceeding 90%, three have hints of PAH emission with AGN 
fractions between 80% and 90%, and the remaining four have 
noticeable contributions from PAH features with AGN frac- 
tions < 80%. For com parison, only 2 of the 13 SMGs from 
iPope et all (l2008bl) and lMurphv et al.1 (l2009h have similarly- 
derived AGN fractions exceeding 60%. One of these 2 
sources is IRBGS from our sample, and the other is an IRBG 
with a 24/im flux density that fell below our flux cut of 700 
^Jy (/24 = 303/xJy). 

As discussed above, the mean offset between the IRAC and 
IRS spectral slopes of the power-law galaxies is only 18%. 
Furthermore, none of the power-law AGN have IRAC and 
IRS slopes that differ by more than ~ 40%, indicating that 
the power-law emission seen at rest-frame wavelengths of 
~ 1 - 3 /im continues up to ~ 12 fim, and possibly beyond, 
assuming a typical redshift of z ~ 2. 

Recall, however, that 3 of the IRBGs in our sample (IRBGS, 
IRBG 10, and IRBGl 1), have far-higher ratios of 24 fim to 8 
/xm flux than the rest (524/^8.0 > 20, see Figure 2), due either 



to extremely high levels of extinction extending into the NIR 
(at z = 2, the 2.7 fim rest-frame emission is redshifted into the 
observed 8 fim band) or to a contribution from the 7.7 fim 
PAH feature. The IRS spectra suggest that both of the above 
scenarios contribute to the enhanced 24 fim flux densities of 
IRBGS and IRBGIO. At redshifts of z = 2.00 and z = 1.92, 
the 7.7 fim PAH feature, clearly visible in each of the spectra, 
serves to enhance the observed MIPS flux. However, star- 
formation is responsible for only ^ 2S% of the 24 fim flux 
density in these two sources. Removal of the star-forming 
flux therefore lowers the observed 24 fim to 8 fim flux ratios, 
but only to values of IS. 9 and 17.6, still well above the ratios 
seen for the remaining sources. The power-law continuum 
at an observed wavelength of 8 fim is therefore likely to be 
suppressed, although shorter-wavelength spectral data would 
be needed to verify this hypothesis. 

The remaining source, IRBGl 1, differs from the others in 
that the AGN accounts for 92% of the MIR flux density, and 
no noticeable contribution from PAH emission is seen in the 
spectrum. In this case, the large 24 fim to 8 fim flux ratio 
most likely results from extreme MIR obscuration. Indeed, 
this source shows one of the largest 9.7 fim opacities in our 
sample, rg 7 = S.O, cor respond i ng to an Ay of 92. S if we as- 
sume the conversion of lDraind (1200 3h . 

5. X-RAY PROPERTIES 

As discussed in the introduction, the study of luminous IR- 
BGs has suffered thus far from a relative lack of sensitive 
X-ray coverage in the fields where these sources have tradi- 
tionally been identified (e.g. the Spitzer First-Look Survey, 
the NOAO Deep Wide Field Survey, and the SWIRE survey). 
As such, the X-ray analysis of IRBGs has lagged behind their 
study in the infrared, and to date, no studies have obtained 
both IRS spectra and deep X-ray imaging of a statistically 
significant sample of IRBGs selected independently of their 
X-ray properties. The deep X-ray coverage and complete IRS 
sampling of our uniformly-selected IRBGs will therefore help 
to bridge the gap between prior X-ray and MIR spectral stud- 
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TABLE 4 
X-RAY Properties 



ID 


Tx 


/{0.5-2 keV) 


7(2-8 keV) 


/(0.5-8 kcV) 


log L2-IO kcV*" 


logA'H'' 


log L2-IO kcV.con'' 




(s) 


(ergs cm~^) 


(ergs s~^ cm~^) 


(ergs s"' cm"'^) 




(ergs s"') 


(ergs s"^) 


IRBGl 








6 29E-15 


43.6 


23 4 


44.3 


IRBG2 


1.57E+06 


2.86E-15 


7.50E-15 


1.06E-14 


44.2 


23-0 


44.6 


IRBG3 


2.23E+05 


<7.07E-17 


1.37E-15 


6.53E-16 


<42.3 


> 23.6 


>43.6 


IRBG4 


1.67E-I-06 


5.58E-15 


7.42E-15 


1.30E-14 


44.2 


22 0+° ' 


44.2 


IRBG5 


1.65E-I-06 


<4.00E-17 


<3.26E-16 


<2.06E-16 


<42.3 




IRBG6 


1.92E-I-06 


2.07E-16 


2.48E-15 


2.52E-15 


43.0 


23 6+0.0 


44.0 


IRBG7 


1.41E-I-06 


1.52E-16 


1.32E-15 


1.40E-15 


42.7 


23 4+^-' 


43.6 


IRBG8 


l.l(,E+01' 


4.28E-15 


1.04E-14 


1.47E-14 


44.3 


22.8 


44.6 


IRBG9 


4.86E+07 = 


<1.86E-16 


<1.10E-15 


5.39E-16'' 


<42.9 




IRBGIO 


5.10E-l-07'= 


1.86E-16 


<4.85E-16 


5.11E-16 


42.9 


< 22.8 


< 43.2 


IRBGl 1 


5.78E-I-07 = 


<6.64E-17 


<8.63E-17 


<1.53E-16 


<42.5 






IRBG12 


7.34E-I-07 = 


<1.74E-16 


<7.76E-16 


<6.49E-16 


<42.3 






IRBG13 


7.47E-I-07 •= 


<1.24E-16 


<7.95E-16 


<4.93E-16 


<43.0 







" Log of the rest-frame, non absorption-corrected 2-10 keV luminosity, calculated using the observed 0.5-2 keV flux, F = 2.0, and an assumed redshift of z = 2.0 
for those sources lacking redshift constraints. Errors are based on the uncertainties in the soft-band flux. 
'° Estimated using the observed flux ratio and an intrinsic photon index of F = 2 
Log of the rest-frame, absorption-corrected 2-10 keV luminosity 

Low-significance X-ray detections: o-(IRBG3, 0.5-8 keV) = 4.8; cr(IRBG3, 2-8 keV) = 5.5; o-(IRBG9, 0.5-8 keV) = 2.1 
° Units of s cm-, nominal exposure time is 200 ks 



ies of these sources. 

5.1. X-ray Detection Fraction 

The X-ray dat a for the CDF-N , CDF - S, E-CDFS-S, an d 
EGS come from Alexander et al. (2003), Luo et al. (2008), 
iLehmeretai. (2005 ). and .Laird et al . (2009), respectively. At 
the position of the 1 3 IRBGs, the X-ray exposm^es range from 
~ 200 ks to ~ 2 Ms, with a mean value of ^ 800 ks. Only 
7 of the 13 IRBGs 50%), however, are formally detected 
in the X-ray (e.g., are detected to high enough significance to 
be included in published X-ray catalogs). Their X-ray fluxes 
in the full (0.5-8 keV), hard (2-8 keV), and soft (0.5-2 keV) 
X-ray bands are given in Table 4. For sources in the EGS, 
we scale the published 0.5-10 keV and 2-10 keV fluxes to 
the bands defined above assuming an observed photon in- 
dex of r = 1.4 (the observed photon index of the typical ob- 
scured AGN, an d thus of the X-ray background, in this en- 
ergy range, e.g., iMarsh all et al... 1980; Gendreau et alj [19951 : 
iHickox & MaSe vitch 200^^ 

For the 6 sources not detected in the catalogs listed above, 
we searche d for faint X-ray co unterparts using the procedure 
outlined in Donlev et alj ( 12005 ). Briefly, we first calculated 
the appropriate 60%, 70%, and 80% encircled energy radii 
(EER) for each source by exposure- weighting t he EER, as 
measured using the lookup tables of iLaird et al.l (i200 9), for 
each individual Chandra observation. After measuring the 
source counts in each of the three apertures, we measured sky 
counts in 10,000 randomly-placed apertures of the same size 
that ( 1 ) did not intersect the source aperture, that (2) lie within 
1' of the source position, and that (3) did not contain a known 
X-ray source. In addition, we corrected the counts in each sky 
aperture to reflect the exposure time measured for the source, 
and only allowed sky apertures in which the exposure differs 
from the source exposure by no more than 15%. We then 
fit a Poisson distribution to the resulting distribution of sky 
counts to determine the significance of any source detection, 
and chose the EER that maximized the source signal. 

Of the 6 formally X-ray non-detected sources, IRBG3 is de- 
tected in both the full and hard bands to 4.8 and 5.5 a above 
the sky, respectively, and 1RBG9 is marginally detected in the 



full band, but only to 2.1a above the sky. The remaining 4 
sources remain undetected in all bands, although 1 (1RBG5) 
lies too close to a known X-ray source to test for very faint 
X-ray emission. We place conservative 3(t upper limits on the 
flux of each X-ray-undetected source by adding any positive 
source counts to a 3(7 upper limit on the measured sky back- 
ground in the 70% EER, assuming the typical observed X-ray 
photon index of an obscured AGN, F = 1 .4. 

5.2. X-ray Obscuration 

For the 8 IRBGs with a formal or weak detection in the 
hard and/or soft X-ray bands, we crudely estimate the intrin- 
sic column density using the hard to soft X-ray flux ratios, 
the measured redshifts (or z = 2 for sources lacking a secure 
redshift measurement), and an assumed intrinsic (e. g., un- 
absorbed) X-ray photon index of F = 2.0 (e.g., ,George et alj 
|2003). This method, commonly used to estimate the column 
densities of faint X-ray sources lacking sufficient counts for 
spectral fitting, is necessarily approximate in nature as it as- 
sumes an X-ray spectrum represented by a single absorbed 
power-law and neglects the potential effects of more com- 
plicated spectral components such as the soft X-ray excess 
and/or Compton reflection. Nonetheless, column densities es- 
timated in this way are broadly consistent with those deter- 
min ed via X-ray spec tral fitting (see, for example. Appendix 
B of iPerola et al.ll2004) . For instance, the column density we 
estimate for IRBG6, log A^H(cm"^)= 23.6, i s in good agree- 
ment with the spectrally-derived value of [Alexander et alj 
(2005), log A^H(cm-2)= 23.8. 

The column density estimates are given in Table 4. To 
place measurement errors on the columns, we take into ac- 
count both the errors in the hard and soft X-ray flux, as well 
as the uncertainty on the redshift. For those sources lacking a 
secure redshift measurement, the redshift error is taken to be 
the standard deviation in the redshift distribution for the full 
sample of IRBGs, CTz = 0.42. Of the 8 sources for which X- 
ray-based estimates could be made, 7 are formally obscured 
with log A^H (cm"^)> 22 and 1 has an upper limit consistent 
with significant obscuration (log A^h (cm"^)< 22.8). 
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From their infrared properties, we know that all of our sam- 
ple contain strong AGN. To place crude column density esti- 
mates on the X-ray-non-detected IRBGs, and to obtain inde- 
pendent estimates for th e X-r ay-detected sources, w e follow 
[Alexander et alj (12008b) and iLanzuisi etal] (12009') and ex- 
plore the observed relationship be tween an AGN's rest-frame 
X-ray and MIR luminosity (e.g., iLutz et all 120041) . The po- 
sition of our IRBGs in MIR/X-ray color space is shown in 
Figure 4, where the rest-frame 6 /im luminosities (i^Lij) have 
been estimated using the best-fit AGN contribution to the IRS 
spectra (which directly sample the rest-frame 6 /im emission 
at z > 1.4) and the rest-frame 2-10 keV luminosities have 
been estimated using the observed soft band flux (0.5-2 keV 
at z = 0, 1.5-6 keV at z = 2) and an assumed observed pho- 
ton index of F = 1 .4, typical of obscured AGN. The shaded 
region in Fi gure 4 gives t he intrinsic observed relation for lo- 
cal AGN (Lutz et al. 2004) and the solid line gives th e intrin - 
sic luminosity-dependent relation of MaioUno et al.l (l2007h . 
who drew their AGN sample from both low (z < 0.2) and 
high (z = 2-3) redshiftfl As shown, the average absorption- 
corrected (e.g. intrinsic) luminosity ratio of the 6 sources 
in our sample with crude X-ray-derived column density es- 
timates is consistent with the Maiolino et al. (2007) relation. 
The ef fect of various o bscuring columns on the observed rela- 
tion of iMaiolino et al.l ( l2007i) is shown by dashed lines, where 
we again assume an intrinsic photon index of F = 2.0. We 
also assume for simplicity that the 6 /im emission from the 
torus is isotropic, i.e. that it does not depend on the ob- 
scuring column or on the angle at which the AGN/torus is 
viewed. If the torus i s optically thick a t MIR wavelengths (e.g . 
iPier & KroUkllT99l iHeckmanI [T99l iBuchanan et al]l2006l) . 
however, an obscured AGN would fall to both lower X-ray 
and lower MIR fluxes, thus increasing somewhat the column 
density estimated via this method. 

For comparison, we also show the absor ption- corrected col- 
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O IR-briqht/opticolly-foint AGN (this work) 



ors of the Type 2 QSOs from 
lel 



jSturm et al.1 (I200M) . all of which 
fall in or above the lLutz et al. 1 (12004 region, and the observed 
colors of the z 2 Compton-thick AGN from lAlexander et"an 
(|2008b), all of which fall below the A^h = lO^'* cm'^ rela- 
tiorQ. Furthermore, we also plot the observed colors of the 
six SMGs with a well-me asured MIR AGN contribution from 
either PopeetaD (l2008bi CI, C3, GN04, GN06, GN07) or 
thi s work (bo t h IRBG 5, which corresponds to source CI 
in iPope et a l.' ('2008b'), and IRBG6 are SMGs). For the 
iPope et aTT i 2008b) sources, we estimate the intrinsic 6/im lu- 
minosity using their Fig. 3 for which they performed spectral 
fitting similar to that d escribed in §4. X-ray fl ux measure- 
ments were taken from [Alexander et al.l (120031) when avail- 
able (CI and GN04), and were calculated as described in §5. 1 
otherwise. Of the 6 SMGs, three (IRBG6, C3, GN04) are 
detected in the X-ray to high significance, two (GN06 and 
GN07) are weakly-detected at 6.2a and 3.3a above the back- 
ground, respectively, and one (C1/IRBG5) remains X-ray un- 



^ We as sume the average IRS spectral slope of our sample, a = —1.94, to 
convert the lMaiolino et al] ilOOlh relation, initially defined at 6.7 fim, to 6.0 
//m. 

' The source SMM J123600+621047 from [Alexander et al .> l2008b) cor- 
responds to source IRBGS in this work. Because of the slightly differ- 
ent procedures used to estimate the rest-frame X-ray and IR luminosities, 
however, our calculated values of log L2-iokev(£rgs s"')= 45.4 ergs s"' 
and log i^L^ (ergs s"')< 42.3 differ slightly from those in [Alexander et al] 
(2p08b): log L2-ioke v(ergs s~*)= 45.3 and log i/L^(ergs s"')< 42.4. For con- 
sistency, we plot the [Alexander et alj fSOOSh) source using the values mea- 
sured here. 
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Fig. 4. — Rest-frame obscured X-ray luminosity vs. rest-frame 6 /xm 
monochromatic luminosity, the values of which are given in Tables 2 and 
4. The IRBGs are shown as circles. The intrinsic relation for local 
AGN I Lutz etal. 2004) and the intrinsic, luminosity-dependent relation of 
[Maiolino et al. 1 ,2007) . which is based on both low and high redshift (z 2) 
AGN, are given by the shaded region and solid fine, respectively. As indicated 
by the red star, the absorption-corrected (e.g. intrinsic) luminosity ratio of the 
six sources in our sample with crude X-ray-derived colum n den sity estimates 
is broadly consistent with the Maiolino et al. 1 200^ and ILutz et al. ( 200^) 
relations (the error bars represent the standard deviation in the distribution 
of X-ray and 6 /^m luminosities). Dashed fines give tfie expected observed 
relations for AGN obscured by column densities of Nyi = 10^^^, 10^"*, and 
2 X 10^* cm"-, calculated assuming an intrinsic X-ray pfioton index of F = 2. 
For comparison, the observed X-ray and intrinsic 6 /im luminosities of the 
Compton-thick AGN from Alexander et al. 1 2008b) are shown as diamonds 
and a sample of SMGs containing AGN drawn both from this work and from 
Pope et al. 1 2008b) are shown as squares. The absorption-corrected values 
for the Type II AGN from Sturm et aL t2006) are shown as triangles. Based 
on this luminosity ratio, 5/13 IRBGs are strong Compton-thick candidates 
with predicted columns of Wh > 2 x 10^"* cm"-, and 5 additional sources lie 
near tfie Compton-tfiick/fieavily obscured boundary (Wh > lO-'' cm"^). 



detected. 

As shown in Figure 4, three of our 13 sources (IRBG2, 
IRBG4, and IRBG8) appear to be obscured yet Compton- 
thin (A^H ^ 10^"' cm"^); the remaining 10 lie near or above 
the Compton-thick boundary (A^e > 10^"* cm"^). This high 
fraction of Compton-thick candidates (~ 80%) i s significantly 
larger than that found by Lanzuisi et al.l ( 120091 18%) thanks 
to the far deeper X-ray data (7x =200 ks to 2 Ms com- 
pared to 7x = 5-70 ks). Furthermore, there is at least rough 
agreement between the column densities estimated above and 
those estimated here. For instance, of the 3 sources with X- 
ray/MIR luminosity ratios indicative of column densities of 
log A^H (cm"^) < 23, all have X-ray-derived column densities 
of 22 < log A^H (cm~^) < 23. Of the remaining three sources 
with X-ray-derived column densities (IRBGl, IRBG6, and 
IRBG7), two are identified as Compton-thick candidates in 
Figure 4 and one (IRBGl) has an estimated column that falls 
just short of log A^H (cm"^) = 24. While their X-ray-derived 
columns (log A^h (cm"^) = 23.4-23.6) place them just outside 
the Compton-thick regime, column densities based on low 
S/N flux ratios are often unable to distinguish Compton-thick 
AGN from their most heavily-obscured, yet Compton-thin, 



IR-Bright/Optically-Faint Galaxies 



9 



counte rparts (see, for example. Figure 3 of [Alexander et al] 
I2008bl) . On the basis of the X-ray and X-ray/MIR-derived 
column densities, we therefore conclude that all of the sources 
in our sample are likely to be X-ray obscured (log A^h (cm"^)> 
22), and as many as 10 (^ 80%) may be Compton-thick. Fur- 
thermore, at least 5 of these sources (^ 40%) lie at predicted 
column densities of > 2 x 10^"* cm~^, and are therefore strong 
Compton-thick candidates. 

One of these strong Compton-thick candidates, 1RBG5, 
was previously identified as a potential Compton-thick quasar 
by both Alexander et al. (2008b) and Georgantopoulos et al. 
U009) and is also one of the SMGs plotted in Figure 4. Its 
heavily obscured nature, however, is not unique among the 
sample of SMGs with measurable AGN contributions in the 
MIR. In fact, 4 of the 6 SMGs plotted in Figure 4 (~ 70%) 
are Compton-thick candidates. Our results are therefore in 
agreement with those of Alexander et al.l (120051) . who found 
that 80% of SMGs containing X-ray-detected AGN are hkely 
to be heavily obscured (log A^h (cm~^) > 23) and potentially 
Compton-thick. This similarity between the column density 
distributions of the IRBG and SMG samples underscores the 
potential relationship between these two classes of objects, 
discussed in more detail in §7 and §8. 

5.3. X-ray Luminosity 

The observed and absorption-corrected X-ray luminosities 
of the IRBGs are given in Table 4. The mean, non-absorption- 
corrected 2-10 keV rest-frame luminosity of the 7 sources 
with soft-band detections (calculated as described in §5.2) is 
log L2-i{)(ergs s~')=43.9, which rises to log L2-io(ergs s"')= 
44.3 when we apply the absorption corrections to the 6 
sources with X-ray-derived A^h measurements. 

As discussed in §5.2, an AGN's intrinsic X-ray luminos- 
ity can also be estimated from its rest-frame 6 fim luminosity 
(see Figure 4). This technique has the added advantage that 
it can be applied to all of the sources in the sample, regard- 
less of whether they are formally detected in the X-ray. As- 
suming a redshift of z = 2 for the 3 sources lacking redshift 
estimates, we find a mean intrinsic 2-10 keV luminosity of 
log L2-io(ergs s~')= 44.6 ±0.4 for the 13 sources in our sam- 
ple. This value is only slightly higher than the estimate made 
above for the sources with X-ray-derived column density esti- 
mates (log L2-io(ergs s"')= 44.3), as would be expected if the 
flux ratio method underestimates the column densities of the 
most heavily obscured sources. Furthermore, Figure 4 sug- 
gests that essentially all of the sources in our sample (with the 
exception of IRBG 12, whose redshift is only z = 1 .04) have in- 
trinsic X-ray luminosities that lie above the canonical division 
between low-luminosity Seyfert galaxies and high-luminosity 
QSOs, log Lx(ergs s"')= 44. The IRBGs in our sample are 
therefore Type 2 QSOs, at least from the X-ray perspective. 

5.4. Correlation between X-ray and IR Properties 

All three of the sources in our sample that are optically- 
thick at 9.7 /im (T9 j^m ^1-5) are Compton-thick candidates 
with MlR/X-ray-estimated column densities of A^h > 2 x 

10^"^ cm"^. The seven additional Compton-thick candidates, 
however, display little to no silicate absorption (rg v^m < 
0.5). While it is possible that the 3 sources lacking solid 
redshift estimates (1 of which is a Compton-thick candidate) 
could have a silicate feature redshifted out of the observable 
band (z > 2.6), the highest confirmed redshift in our sample 



is that of IRBG13, z = 2.6, a source with a prominent absorp- 
tion feature. Given the average sample redshift of z = 2 and 
the low frequency of silicate absorption among the sources in 
our sample with confirmed redshifts, it therefore seems un- 
likely that the three sources lacking redshift information (and 
observable silicate features) lie at z > 2.6 and have significant 
silicate absorption. 

It therefore appears clear that while silicate absorption is 
more likely amongst heavily X-ray-obscured IRBGs, strong 
silicate absorption need not accompany X-ray obscuration. 
A similar conclusion was drawn by Brand_et al. ( 2008) for a 
sample of luminous (/24 > 2 r nJy) IRBGs in the N OAO Deep 
Wide-Field Survey, as well as ISturm etaTI (l2005h . who show 
that bright (Fis^m > 0.3 mJy) X-ray-selected Type 2 AGN 
can show unabsorbed power-law continua in the IR despite 
heavy obscuration in the X-ray. This apparent discrepancy be- 
tween X-ray and MIR measures of obscuration could be due 
to an additional component of absorbing material t hat blocks 
the X -ray, but not the IR, emission region (e.g., IShi et al.l 
I2OO6I) . Alternatively, it could result from the filling-in of the 
silicate feature by an additional component of extended hot 
dust that lies beyond the obscuring torus, arising, for exam- 
ple, from the ionization cones as in NGC 1068, the prototyp- 
ical Compt on-thick AGN in the local Universe (iMason et al.l 
12006- .Efstathiou et al.|[T995h . 

6. STAR-FORMATION RATES 

For the 7 IRBGs with detectable PAH emission (i.e. AGN 
fractions < 90%), we place limits on the star-formation 
rates (SFRs) using (1) the star-forming contribution to the 
observed-frame 24 emission (in lieu of direct estimates 
based on the weakly-detected PAH emission features) and (2) 
the observed radio emission. 

6.1. 24 emission 

To measure the 24 /xm flux density arising from star- 
formation alone, we convolve the star-formation component 
of our IRS fit with the MIPS 24 passband. We find that in 
the sources with detectable PAH emission, ^ 25% of the to- 
tal 24 /im flux density can be attributed to star-fo rmation. We 
then use the redshift-dependent scaling relation ofl Rieke et al.l 
(2009) to convert the observed MIR flux to a SFR, assuming 
the median sample redshift of z = 2 for those sources lacking 
secure redshift estimates (for further details, see Rieke et ai] 
|200^. With the exceptions of IRBG 13, whose anomalously 
high SFR estimate of 2 x 10^ Mq yr"' will be discussed be- 
low, and IRBG12, whose estimated SFR is ~ I5OM0 yr"', 
the resulting SFRs range from - 1000-2000 M© yr"'. No 
meaningful constraints can be placed on the SFRs of those 
sources with AGN contributions in excess of 90% (SFR< 
7000-5 X 10"* Mq yr-'). 

These estimates, however, are based on a relationship de- 
rived from local starburst, LIRG, and ULIRG templates. At 
high-redshift, luminous star-forming galaxies exhibit strong 
PAH features typical of less-luminous local galaxies, due per- 
haps to extended star-formation or lower metallicity (e.g., 
Saiinaetal. 2007; Papovich et al. 2007; Pope et al. 2008b[ 
Rigbvetal. 2008; FaiTahetal. 2008; Mur phy et al.l i2009t 
Menendez-Delmestre et al. 2009). Local templates are there- 
fore likely to overestimate the intrinsic 24 /im flux (and there- 
fore the SFR) of high-z ULIRGS significantly at L{TIR) = 
10'^ Lf7^, the typ i cal luminosities o f luminous IRBGs (see 
iRieke et al.1 120091: iTvler et all 120091) . Additional uncertain- 
ties in the derived SFRs arise from the lack of luminous 
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L(TIR) > 2 X 10'2 Lq ULIRGS in the local Universe from 
which to calibrate the SFR relation. The IR-derived SFRs of 
the IRBGs should therefore be taken as rough upper limits on 
the true SFRs. 



6.2. Radio emission 

To place an independent constraint on the SFRs, we there- 
fore consider the 1.4 GHz radio flux density. Of the 7 IR- 
BGs with detectable PAH_emission, 6 have radio counter - 
pa rts in the catalogs ofjRichards (2000'), llvison et al] (l2007h . 
or lKellermann et all (I2008I) . Assuming a radio spectral index 
of Of = -0.7, where f^, (x ly", we K-correct the o bserved ra- 
dio flux and use the relation of lRieke et al.l (120091) to estimate 
the SFR. With the exception of IRBG13, whose radio-derived 
SFR falls just short of the IR-derived SFR, the radio SFRs ex- 
ceed the IR SFRs by factors of ^ 1.2-10, suggesting that a 
significant fraction of the radio flux density comes not from 
star-formation, but from AGN activity. Therefore, the radio- 
derived SFRs again represent only an upper limit on the true 
star-formation activity. 

To isolate the star-forming contribution to the radio emis- 
sion, we therefore turn to the radio-infrared correlation, char- 
acterized by the parameter ^24 = log(/24/im//i.4GHz)- By us- 
ing this correlation together with the measured star-forming 
contribution to the 24 /im flux density to calculate the star- 
forming radio flux density and then the SFR, we again uti- 
lize the MIR emission as a proxy for the SFR. This method, 
however, allows us to calibrate the SFR using the observed 
radio/SFR correlation in place of the 24/xm/SFR correlation 
discussed above. 

A number of different studies have attempted to quan- 
tify (7?4 and have returned values ranging from 0.52 
dBeswicket al.l l2008) to 1.39 (Boyle etal. 2007). Here, we 
adopt the luminosity-dependent ^24 relation of Rie ke et al] 
(|2009) which is based on the local IRAS bright galaxy sam- 
ple (BGS) and for which ^24 = (-1.275 ± 0.756)+ (0.224 ± 
0.066) X log L(TIR)/Lq at log L(TIR)/Lq > 11. At 
log L{TIR)/Lq = 1 3.0, the typical total infrared luminosity 
of luminous IRBGs (iTvler et al.l l2009), ^24 = 1-64. 

The observed value of q24 also varies with redshift as the 
IR spectral features pass through the 24 /im band. To cor- 
rect for this variation, we a dopt the observed ^24 redshift 
evolution of llbar et aLl (120081) . who found that ^24 = (0.85 ± 
0.01) + (-0.20±0.01) xz. Not surprisingly, this observed red- 
shift evoluti on is reproduced b y the log L{TIR)lLp i ~ 12.0 
templates of lRieke et al.l (12009 ). After scaling the Ibar et al. 
(I2OO81) relationship upwards to match the .Rieke et al. (2009.) 
result at z = 0, we used the resulting correlation and our esti- 
mate of the star-forming 24 /im flux density to determine the 
expected radio emission from star-formation. 

We find that star-formation accounts for 3% - 23% of the 
total radio emission in our sources. Again with the excep- 
tions of IRBG12 and IRBG13, whose estimated SFRs are - 
1600 Mq yr"' and 20 Mq yr"', respectively, we derive SFRs 
of ^ 200-400 M© yr"' using the star-formation-dominated 
radio emission. Because the high-z ULIRGS/HyperLIRGS in 
our sample are likely to have star-forming emission similar 
to that of local LIRGS/ULIRGS (see §6.1), however, the true 
intrinsic value of ^24 is likely to be lower than the value as- 
sumed above. The estimated SFRs can be therefore be taken 
as lower-limits, subject of course to the intrinsic scatter in the 
radio/IR correlation. 



6.3. The high SFRs of luminous IR-bright/optically-faint 
galaxies 

To summarize, 5 of the 7 IRBGs in our sample with de- 
tectable PAH emission and/or radio counterparts have esti- 
mated SFRs in the range of - 200 to 2000Mq yr"'. The 
SFR of IRBG12 (- 20-150 Mq yr"') is lower than the oth- 
ers because of its significantly lower redshift (z = 1.04). The 
SFR of IRBG13, however, is anomalously high: 1.6-2.0 x 
10'* Mq yr~\ This overestimation likely stems from the fact 
that both an AGN+SF fit and an AGN-only fit reasonably re- 
produce the IRS spectrum (x^ = 0.79 and 0.85, respectively). 
While potential AGN-only fits likewise exist for many of the 
sources in our sample, IRBG13 has the smallest offset in the 
reduced values of the AGN+SF and AGN-only fits and is 
therefore the source for which the origin of the IR emission is 
most uncertain. 

The SFRs derived above are clearly approximate in na- 
ture. Nevertheless, two of the 7 sources examined above 
(IRBG5 and IRBG6) are known SMGs, heavily dust- 
enshrouded sources wi t h typical SFRs of 500 - 2000Mq yr"' 
("Alexander elal] 120051; ISwinbank efa[\ 120041) . As their IR 
properties are typical of the rest of our sample, it therefore 
seems likely that our estimated SFRs of 200-2000Mo yr"' 
reasonably approximate the high level of star formation 
present in luminous IRBGs. Further more, our est i mated SFRs 
are in agreement with the results of ITvler et al.l (120091) . who 
find via far-infrared measurements that approximately a third 
of IRBGs similar to those selected here have indications of 
star formation at a level that produces a substantial fraction 
of their bolometri c luminosities, which are typically 10'^ Lq 
( ITvler et al.112009 '). Finall y, such SF Rs are also in agreement 
with the predictions of Nar ayanan e t al. ( 2009a), who find that 
luminous IRBGs are best modeled by gas rich mergers un- 
dergoing rapid AGN growth and/or star-formation at rates of 
-500-lOOOMoyr-i. 

Presumably, the objects with detected star formation repre- 
sent the high end of a continuous luminosity function, where 
the lower levels of star formation are masked by the AGN 
outputs. That is, it is likely that these sources are in general 
characterized by star formation in the high-LIRG to ULIRG 
range, even though their MIR (14-37 /im) luminosities are 
AGN-dominated. 

7. HOST GALAXY MORPHOLOGIES 

Because of the ^ 30% overlap of IRBG and SMG samples, 
as well as their similar redshift distributions, surface densi- 
ties, and clustering properties, several authors have suggested 
that luminous IRBGs may evolve from SMGs (e.g..'De y et al. 
2008; Pope et al. 2008a; Brodwin et al. 2008; Bussman n et al. 
2009^ JVaravanan et al. 2009b a). Furthermore, early stud- 
ies of SMGs suggested that 61% ±21% may be in the early 
stages of a major merger (iConsehce et al.ll2003bi) . although, 
on the basis of asymmetry measurements alone, more recent 
work has suggested that SMGs are no more likely to be un- 
dergoing major mergers than lo wer-luminosity high-redshift 
galaxies (ISwinbank et al.ll2010l) . While this work has called 
into question the uniqueness of the role that mergers play in 
these systems, the merger-driven origin of both SMGs and lu- 
minous IRBGs i s sup ported by the GADGET2 simulations of 
iNaravanan et al.l (l2009hia) . who find that both samples natu- 
rally arise from gas-rich mergers, with SMGs representing the 
early, star-formation dominated phase of the merger, and IR- 
BGs representing the later phase of final coalescence in which 
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Fig. 5.^ goods and AEGIS ACS ('-band images of the 8 IRBGs with deep ACS coverage (top. image diameter of 2") and the 8 GOODS-N SMGs from the 
literature with similar redshifts and r'-band magnitudes (bottom, image diameter of 4"). The redshift and IRS-derived AGN fraction are given in the top left and 
light comers of the image, and the asymmetry and concentration indices for each galaxy are given in the bottom left and right corners of each image, respectively. 
The redshifts and AGN fractions for the SMGs come from Pop e et al. 1 2006. .20(j8b.) and, Murphy et al„ (2009). The circle gives the measured Petrosian radius, 
and the cross-hairs indicate the rotation center chosen by the asymmetry algorithm. Because it extends beyond the edge of the image, we are unable to estimate 
a reliable Petrosian radius or morphological indices for IRBG7. 



Star-formation and obscured bl ack hole growth both proceed 
at rapid rates (see Figure 1 in Ho pkins etai]|2008h . If this 
scenario is correct, luminous IRBGs should therefore exhibit 
signs of recent merger activity, although one might expect 
them to be more dynamically relaxed than their SMG pre- 
decessors. 

7.1. Comparison between IRBGs and SMGs 

Previous morphological studies of IRBGs have used tar- 
geted HST and Keck AO observations in the Bootes Field of 
the NOAO Deep Wide-Field Survey to (1) measure their typ- 
ical sizes and shapes via Sersic fitting and GALFIT galaxy 
decomposition and to (2) characterize their morphologies 
on the basis of the Gini, M20, and concen tration indices 
jBussmann et alj 120091: [Melbourne et all l2009i) . They con- 
clude that IRBGs are preferentially found in disk-like galax- 
ies, and only '^15% show multiple resolved components in- 
dicative of ongoing mergers, far lower than the ^ 50% in- 
teraction rate found for the average z = 2 ULIRG and the ~ 
60% merger rate estimated for SMGs dConseli ce et al.ll 2003bl: 
Dasvra et all 120081: iBussmann etal] 120091: [Melbourn e et alj 



2009ir However, IRBGs do appear to have smaller pro 
jected sizes (rpetrosian = 0.5"- 1.5") than SMGs (rpen-osian = 



0.5"-2.5"), suggesting that they may in fact r e present a later , 
more relaxed phase of a major merger ( Bussmann et al. 2009). 

To expand upon these studies and further test whether IR- 
BGs are the product of major-mergers, we directly com- 
pare below the morphologies of IRBGs and SMGs in a self- 
consistent manner Because IRBG5 and IRBG6 are confirmed 
SMGs, we treat these two sources separately in the analysis 
below. We also note that because the EGS does not yet con- 
tain deep sub-mm coverage, it is possib le that one or mor e 
of the EGS IRBGs could also be SMGs. iPope et alj (l2008al) . 
however, found an SMG fraction of only 30% amongst lu- 
minous IRBGs. As the fraction of SMGs in our sample of 8 
IRBGs with deep ACS coverage is akeady ^ 25%, we there- 
fore expect little to no additional contamination by SMGs in 
our morphological sample. 

We select as a compari son sample of S MGs those sources 
in the GOODS-N field (iPope et alj|2005l) with redshifts and 
apparent magnitudes consistent with those of the IRBGs: 
1 < z < 3 and / < 26, and confirm that none are also IR- 
BGs (although GN25 has a high /24 //r of - 900: lCowie etaP 
[2004; Pope et al. 2006). The observed-frame /-band images 
of the 8 iRBGs and 8 SMGs with deep ACS imaging from 
GOODS or the EGS are shown in Figure 5. 



12 



DONLEY ET AL 



7.1.1. Merger Fraction 

To quantitatively compare the morphologies of these two 
samples, we calculate the Petrosian radii and concentra- 
tion (C ) and as ymmet ry (A) indices using the approach of 
Conseli ce et al. (120001) and the error-estimation method of 
Shi et al.i (i2009i) . and confirm excell ent agre ement between 
our measured parameters and those of I Shi et al. (2009) for the 
two sources that lie in both samples: IRBG5 and 1RBG6 (Y. 
Shi 2009, private communication). 

Of these parameters, the asymmetry index is the most rel- 
evant as it tends to be high (A > 0.35 in the local Universe) 
when merging galaxies ar e undergoing either the i r first pass 
or their final coalescence (IConselice et alj|2003at iLotz et al] 
I2OO8). We caution, however, that this index is designed to 
measure the asymmetry within the Petrosian radius (rpeu-osian) 
of the primary source. If the merging galaxies are separated 
by a large distance (as occurs, for example, during the merger 
stages between first pass and final coalescence), the secondary 
source will fall outside of the Petrosian radius of the primary 
source and the index will be low. Furthermore, surface bright- 
ness dimming causes the observed asymmetry to decrease 
with increasing redsh ift. While the co rrection for this effect is 
somewhat uncertain, IConselice et al.l ( I2b05i) estimate an off- 
set of AA = -0.15 between z = Q and z = 2 for irregular galax- 
ies. As the ShietaLl ("2009^ error method results in measured 
asymmetries 0. 05 lower than those fou nd by the minimum er- 
ror method of IConselice et alJ (120001) . we therefore adopt a 
merger index of A > 0. 15 for our z ~ 2 sampl43 

The measured concentration and asymmetry indices are 
shown in Figure 6. By these measures alone, the IRBGs and 
SMGs appear to have similar morphologies. Furthermore, all 
of the sources with clear double nuclei have A > 0.15, as ex- 
pected. Several SMGs, however, do indeed have nearby coun- 
terparts that lie outside the Petrosian radius of the primary 
source. If we therefore count as mergers all sources with ei- 
ther A > 0.15 or nearby counterparts that fall outside rpetrosian 
(e.g. GN22,GN20.2,GN20), we find that 5 of the 7 (- 70%) 
IRBGs and 5 of the 8 (-- 60%) SMGs are merger candidates. 

It is also worth noting that two of the three IRBGs which are 
optically thick at 9.7 ^m (IRBGll and 1RBG13) are strong 
merger candidates with large asymmetries and clearly dis- 
turbed morphologies. The third, 1RBG7, falls too near the 
edge of the ACS image to accurately calculate the asymme- 
try. Nonetheless, it appears to lie at the edge of a large patch 
of diffuse emission, potentially due to an ongoing merger (see 
Figure 5). As discussed in §5.4, these three optically-thick 
sources are the only Compton-thick candidates with strong 
silicate absorption. It therefore appears plausible that ongoing 
mergers are at least partially responsible for the high degree 
of MIR obscuration in these heavily X-r ay-obsc ured IRBGs. 
In their sample of local ULIRGs, Veill eux et alJ (|2009) simi- 
larly find that the depth of the 9.7 /im silicate feature is gener- 
ally largest during the close pre-merger and merger stages. 
They caution, however, that this correlation suffers from a 
large degree of scatter, likely due to variations in both the ini- 
tial conditions (e.g., the structure, surface mass density, and 
gas fractions of the infalling galaxies) and specific character- 
istics (e.g., orbital geometry and mass ratio) of the mergers 
themselves. This scatter may explain why the three remain- 

We note that'Shi et al.' ("2009) find a slope (AA/Az) between z, = and 
z = I that is twice that of Conselice et al. 1 2005) when using local LIRGs as 
input. It therefore remains possible that galaxies in our sample with A < 0.15 
are undergoing mergers as well. 
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Fig. 6. — Asymmetry vs. concentration indices for the IRBGs (black 
squares), SMGs (red triangles), and sources that meet both criteria (green 
stars). The mean properties and standard deviations of the IRBG and SMG 
samples are shown by large symbols. As noted in §7.1.1, the measured asym- 
metry will be low when merging galaxies lie at large separations, as is the 
case for many of the SMGs. To calculate a merger fraction, we therefore con- 
sider both the measured asymmetry and the fraction of sources with nearby 
neighbors. 

ing merger candidates (IRBG5, IRBG6, and IRBGIO) show 
little to no MIR obscuration (rgj^im = 0.0-0.3). 

7.1.2. Physical Size 

While a majority of both SMGs and IRBGs appear to be 
sufficiently disturbed to be classified as major mergers, these 
two samples do appear to differ i n their physical size, as in- 
ferred by B ussmann et al.l (120091) . A histogram of the Pet- 
rosian radii in physical units of kpc is shown in Figure 7. 
The mean rpetiosian of the SMGs, 6.8 ±3.7 kpc, is nearly twice 
that of the IRBGs, 3.7 ± 1 .8 kpc. If we exclude the 3 lowest- 
redshift SMGs with z < 1.4, the mean rpetrcsian of the SMGs 
drops to 4.5 ± 2.0 kpc, suggesting that the apparent offset in 
physical size may b e driven at least in part by evolution in 
the SMG population. ISwinbank et al.l (I201Q) . however, detect 
no significant size evolution in either the optical or NIR radii 
of a sample of 25 SMGs. Furthermore, they find a mean op- 
tical (NIR) Petrosian radius of 6.9 ± 0.7 kpc (7.7 ± 0.6 kpc) 
for their SMG sample, consistent with our results and again a 
factor of ^2 higher than observed for the IRBGs. 

This apparent offset in size may be a natural byproduct of 
the merger hypothesis, as the Petrosian radius of a merger 
remnant peaks strongly during both the first pass and the fi- 
nal merger stages and can easily obtain a value twice that of 
the pre-merger, maximal separation (when only one galaxy is 
contained within the Petrosian radius), and post-merger stages 
dLotz et alj|2008i J. Lotz 2009, private communication). An 
offset in rpetrosian of a factor of < 2 is therefore consistent with 
the hypothesis that SMGs are observed in the first stages of a 
major merger whereas IRBGs represent a post-merger phase 
following the final coalescence. What remains more difficult 
to explain, however, are the simultaneously large asymme- 
tries and small radii of the IRBGs, as these two parameters 
should both peak at approximately the same stages during the 
merger. Larger samples of such sources are clearly needed to 
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Fig. 7. — Petrosian radii of the IRBGs (green histogram, angle=-45°), 
SMGs (red histogram, aiigle=45°), and sources that meet both criteria (blue 
histogram, angle=0°). The mean radii (and standard deviations) for the three 
samples are shown at the top of the plot. 

resolve this apparent discrepancy. 

8. DISCUSSION 

Several authors have suggested that IRBGs and SMGs 
may be evolutionarily related, each representing a different 
phase in the major-merger scenario originally proposed by 
jX988) (e.g. Dev et al. 2008; Pope et al. 2008a 



Sanders et al 
Copi)in et al.l 



2010t iNarayana n et al.ii2009biat iBrodwin et al.l 



2008h . To test this hypothesis, we explore below the agree 



ment between the properties of the IRBG and SMG samples, 
as measured here and in the literature, and the major-merger 
model. 

The major-merger (e.g., 3:1 or greater mass ratio) scenario, 
as presented by Hopkins et al. (2008), proceeds as follows. 
During the merging galaxys' first pass, star-formation activ- 
ity begins to increase with concurrent AGN activity occur- 
ring only under certain orbital geometries and/or when the 
host galaxys' discs are particularly unstable. As the galaxies 
undergo final coalescence, however, massive inflows power 
both strong starburst activity and heavily-obscured, yet rapid, 
black hole growth. Eventually, feedback from the black hole 
and from supernovae disperse the remaining gas in a brief 
"blowout" phase characterized by a dust-obscured, yet Type 
1, AGN. Finally, when the gas and dust have been fully dis- 
persed, the source appears as an unobscured quasar (unless 
viewed through the canonical torus). The general picture 
is therefore one in which a major-merger fuels both star- 
formation and black hole growth (thus leading to the ob- 
served Mbh/c relation), with star-formation dominating at 
early times, and AGN activity at later times. 

Before continuing, it is worth noting that this model is 
almost cert ainly an oversimplification. As emphasized by 
IVeilleux eT al. (2009), the internal structure of the merging 
galaxies, the orbital parameters (e.g., prograde or retrograde), 
the mass ratio, and sub-resolution physics will ultimately de- 
termine the evolutionary path followed by any one merg- 
ing galaxy pair, leading to variations on the model presented 
above and scatter in the observed correlations. Nonetheless, 



the general picture presented above should hold on average if 
this process is responsible for the SMG and IRBG phenom- 
ena. 

In the context of this model, SMGs are assumed to rep- 
resent the early, star-formation dominated phase of the ma- 
jor merger, and IRBGs the later, AGN-dominated phase. To 
first order, this progression is supported by the fact that while 
~ 85% of S MGs ar e star-formation dominated in th e MIR 
dPope et al.l l2008b: Hainline et al.ll2009l: ICoppin et aLl lIOlO). 
and nearly all luminous IRBGs are AGN-dominated in the 
MIR, there is an 30% overlap in the SMG and IRBG popu- 
lations, at least for samples extending down to f24 > 100 /iJy 
dPope et al.ll2008ah . Furthermore, the SMGs that also meet 
the IRBG criterion tend have the largest AGN contributions, 
and are therefore Ukely candidates to be undergoing a transi- 
tion from star-formation to AGN dominance (e.g.. Pope et al.l 
l2008aiCoppin et al.ll20T(ih . 

In addition, while SMGs tend to be star-formation- 
dominated in the M IR, --28-50% host obscured AGN 
([Alexander et al.l I2005L also see Figure 4) as predicted for 
the subset of mergers with initial conditions and orbits con- 
ducive to early (obscured) black hole growth. Furthermore, at 
least ^30-4 0% of luminous IRBGs have ULIRG-level SFRs 
(iTyler et al.i 2009, this work). As rapid star-formation is pre- 
dicted throughout the obscured AGN phase of a major merger, 
this too is consistent with the major merger scenario. 

T his general picture is also supported by th e work 
of [Alex ander et al. (2008a), Dey elaH ([20081). and 
'Brodw in et ai.i (120081) . For instance, [Alexander et al.l 
(2008 al) find that black hole growth in SMGs lags behind 
the growth of the host galaxy, contrary to what is seen in 
high-redshift optical (e.g., unobscured) and radio-selected 
samples of quasars (e.g., [Peng et all [2006t [McLure et al.l 
I2OO6). This finding suggests that SMGs are at an earlier 
evolutionary stage than their unobscured, AGN-dominated 
counterparts, ag ain in agreement with the merger scenario. 
More generally, iDev et all ([20081) find that the surface den- 
sities and redshift distributions of /24 > 300 /zJy IRBGs (at 
least ^ 60% of which are likely to be AGN-dominated in the 
MIR,|Pope et al. 20083) a re comparable t o those of luminous 
SMGs (igso^m > 6 mJy, Chapman et al. 2005; Coppi n et aP 
2006). Similarly. , Brodwin et al..(2008) analyze the clustering 
properties of IRBGs and find that /24 > 300 /xJy I RBGs have 
a corr elation length comparable to that of SMGs dBlain et all 
[2004 . The most luminous (/24 > 700 ^Jy) IRBGs, however, 
may be more strongly clustered, and thus lie in richer envi- 
ronments, than both their lower MIR-luminosity counterparts 
and the typical SMG. 

While no one of these arguments confirms the poten- 
tial merger origin for SMGs or IRBGs, collectively, the 
results suggest that the merger scenario first proposed by 
Sanders et al. (1988) may provide an acceptable explanation 
for the properties of these two samples. Ideally, one might 
hope to either confirm or reject this hypothesis on the basis 
of the observed morphologies. Unfortunately, however, the 
morphological properties of these samples remain amongst 
their most uncertain characteristics. As discussed in §7.1, the 
asymmetry measure, commonly used as a merger indicator in 
the local Universe, suffers at high redshift from surface bright- 
ness dimming, a (1 + z)^ effect, and is insensitive to interme- 
diate merger stages in which the gala xies lie at large s epara- 
tions. When using this measure alone. ISwinbank et al.[ ([2010[) 
therefore conclude that SMGs are no more likely to appear as 
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mergers than a bolometrically less luminous sample of high 
redshift star-forming galaxies. Optical (e.g., /-band) mea- 
surements of this and other objective morphology indicators 
are also complicated by the fact that at z ^ 2, the observed- 
frame optical emission probes the the rest-frame UV emis- 
sion, which tends to be more strongly clumped. However, 
when morphologies of AGN-dominated IRBGs are measured 
at longer wavelengths, the contribution from the unresolved 
AGN e mission, now less obscured, can likewise b ias the re- 
sults (Buss mann et al.ll2009l: rMelbourne et al. P2009). 

Prior studies of IRBG morphology have therefore turned 
to the fraction of multiple resolved components as a measure 
of ongoing mergers, and have fou n d a rate of only ~ 15% 
for IRBGs ( iBussmann et al.1 120091: [Melbourne etal] 120091) . 
far lower than is seen in local and high redshift samples of 
ULIRGs. However, this measure too is biased, and will not 
identify sources in the latest stages of a merger We therefore 
chose to combine the asymmetry and near-neighbor/multiple- 
component methods, and conclude that the majority of both 
SMGs and IRBGs are likely to be undergoing mergers. Fur- 
thermore, th e smaller physical size of the IRBGs, a trend 
first noted by Bussm ann et all (12009 ) and confirmed in §7. 1 .2, 
suggests that they are likely in a later, more relaxed merger 
stage, a hypothesis also supported by the AGN and star- 
forming characteristics presented above. Clearly, however, 
one would like to confirm this finding using a far larger, 
deeper, and higher-resolution sample, one that will be pro- 
vided by the upcoming CANCELS (Cosmic Assembly Near- 
IR Deep Extragalactic Legacy Survey) Hubble WFC3 NIR 
imaging survey of the GOODS, COSMOS, EGS, and UDS 
fields. 

9. SUMMARY 

In summary, we present the X-ray, star-forming, and mor- 
phological properties of a sample of luminous (/24 > 700 /iJy) 
IR-bright/optically-faint galaxies (IRBGs, /24//R > 1000) se- 
lected in deep X-ray fields. Our major findings are as follows: 

• The infrared colors and IRS spectra confirm that all of 
the /24 > 700//Jy IRBGs are AGN-dominated. 

• All of the sources appear to be X-ray obscured, and 
4 remain X-ray non-detected despite deep X-ray cov- 
erage. As many as ^ 40% are strong Compton- 
thick candidates, and ^ 40% more appear to lie near 
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the Compton-thick/heavily obscured (Nh = 10^^ cm ^) 
boundary. 

• The mean intrinsic X-ray luminosity of the sample is 
log Lo.5-8(ergs s"')~ 44.6 ±0.4, and all but one source 
has a MlR-derived intrinsic X-ray luminosity in excess 
of log Lo,5-8(ergs s"')= 44, the canonical dividing line 
between Seyfert galaxies and QSOs. These sources are 
therefore Type-2 QSOs from the X-ray perspective. 

• In general, the sources with the largest X-ray obscura- 
tion are those most likely to exhibit strong silicate ab- 
sorption. Furthermore, at least two of the three IRBGs 
that are optically thick at 9.7 /^m are likely to be under- 
going mergers. However, silicate absorption need not 
accompany X-ray obscuration, and merger candidates 
are observed that lack strong silicate features. 

• 30-40% of IRBGs have star formation rates in the 
ULIRG range: presumably the rest of the sample also 
has relatively high levels of star formation that are 
masked by the AGN outputs. 

• As many as ^ 70% of the IRBGs in our sample show 
asymmetries indicative of mergers, similar to the inci- 
dence of mergers in SMGs at similar redshift. However, 
the IR-bright/optically faint sources tend to be more 
compact than the SMGs, suggesting that they may be 
in a later stage of merging. 

• These characteristics are consistent with the proposal 
that these objects represent a later evolutionary stage 
following soon after the star-formation-dominated one 
represented by the SMGs. 
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